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1. Introduction 

The study and analysis of fluctuations are an essential method to character- 
ize a physical system. In general, one can distinguish between several classes 
of fluctuations. On the most fundamental level there are quantum fluctu- 
ations, which arise if the specific observable does not commute with the 
Hamiltonian of the system under consideration. These fluctuations proba- 
bly play less a role for the physics of heavy ion collisions. Second, there are 
"dynamical" fluctuations reflecting the dynamics and responses of the sys- 
tem. They help to characterize the properties of the bulk (semi-classical) 
description of the system. Examples are density fluctuations, which are 
controlled by the compressibility of the system. Finally, there are "triv- 
ial" fluctuations induced by the measurement process itself, such as finite 
number statistics etc. These need to be understood, controlled and sub- 
tracted in order to access the dynamical fluctuations which tell as about 
the properties of the system. 

Fluctuations are also closely related to phase transitions. The well 
known phenomenon of critical opalescence is a result of fluctuations at all 
length scales due to a second order phase transition. First order transitions, 
on the other hand, give rise to bubble formation, i.e. density fluctuations 
at the extreme. Considering the richness of the QCD phase-diagram as 
sketched in Fig^ the study of fluctuations in heavy ions physics should 
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Fig. 1. Schematic picture of the QCD phase-diagram 



lead to a rich set of phenomena. 

The most efficient way to address fluctuations of a system created in 
a heavy ion collision is via the study of event-by-event (E-by-E) fluctua- 
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tions, where a given observable is measured on an event-by-event basis and 
the fluctuations are studied over the ensemble of the events. In most cases 
(namely when the fluctuations are Gaussian) this analysis is equivalent to 
the m eas urement of two particle correlations over the same region of accep- 
tance^! Consequently, fluctuations tell us about the 2-point functions of 
the system, which in turn determine the response of the system to external 
perturbations. 

In the framework of statistical physics, which appears to describe the 
bulk properties of heavy ion collisions up to RHIC energies, fluctuations 
measure the susceptibilities of the system. These susceptibilities also de- 
termine the response of the system to external forces. For example, by 
measuring fluctuations of the net electric charge in a given rapidity inter- 
val, one obtains information on how this (sub)system would respond to 
applying an external (static) electric field. In other words, by measuring 
fluctuations one gains access to the same fundamental properties of the 
system as "table top" experiments dealing with macroscopic probes. In the 
later case, of course, fluctuation measurements would be impossible. 

In addition, the study of fluctuations may reveal information beyond 
its thermodynamic properties. If the system expands, fluctuations may be 
frozen in early and thus tell us about the properties of the system prior 
to its thermal freeze out. A well known example is the fluctuations in the 
cosmic microwave background radiation, as first observed by COBE^. 

The field of event-by-event fluctuations is relatively new to heavy ion 
physics and ideas and approaches are just being developed. So far, most of 
the analysis has concentrated on transverse momentum and the net charge 
fluctuations. 

Transverse momentum fluctuations should be sensitive to tempera- 
ture/energy fluctuations | 57 | 53 | These in turn provide a measure of the 
heat capacity of the system E31 Since the QCD phase transition is asso- 
ciated with a maximum of the specific heat, the temperature fluctuations 
should exhibit a minimum in the excitation function. It has also been argued 
| 55 | 56 | these fluctuations may provide a signal for the long range fluc- 
tuations associated with the tri-critical point of the QCD phase diagram. 
In the vicinity of the critical point the transverse momentum fluctuations 
should increase, leading to a maximum of the fluctuations in the excitation 
function. 

Charge fluctuations EES on the other hand, are sensitive to the frac- 
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tional charges carried by the quarks. Therefore, if an equilibrated partonic 
phase has been reached in these collisions, the charge fluctuations per en- 
tropy would be about a factor of 2 to 3 smaller than in an hadronic scenario. 

In this review, we will systematically examine the principles and the 
practices of fluctuations such as the momentum and the charge fluctua- 
tions as applied to the heavy ion collisions. In doing so, various concepts of 
"average" need to be introduced. They are: 

(i) Thermal average 

(ii) Event-by-event average as a whole 

(iii) Averages over various parts of event-by-event average such as clus- 
ters and particles emitted by clusters. 

When necessary, we will use subscripts to distinguish these various averages. 
However whenever it is clear from the context, or if a relation is of a general 
nature such as the definition of a variance, we will simply use the symbol 
(...) to denote the average. 

The rest of this review is organized as follows. In section 2, we briefly 
review the basic concepts of thermal fluctuations and its relevance to heavy 
ion physics. In section 3, we examine possible sources of fluctuations other 
than the thermal ones and their effect on experimental measurements. In 
section 4, the relationships between the underlying correlation functions 
and charge fluctuations and also balance functions are established and 
their physical meaning made clear. In section 5, various ways of measur- 
ing fluctuations proposed so far in literature are briefly reviewed and their 
inter-relationships established. In section 6, past and current experimental 
results are examined in the light of the preceding discussions. We conclude 
in section 7. 

We regret that due to the lack of space and also expertise on the authors' 
side, such interesting topics as Disoriented Chiral Condensate or Hanbury- 
Brown-Twiss effects couldn't be discussed in detail in this review. Also for 
the same reason, it was impossible for us to cover all the relevant refer- 
ences although we tried our best. If there is any such glaring omission, we 
apologies to the authors. 
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2. Fluctuations in a thermal system 

To a good approximation the system produced in a high energy heavy ion 
collision can be considered to be close to thermal equilibrium. Therefore let 
us first review the properties of fluctuations in a thermal system. Most of 
this ca n be found in standard textbooks on statistical physics such as e.g. 
RefP and we will only present the essential points here. 

Typically one considers a thermal system in the grand-canonical en- 
semble. This is the most relevant description for heavy ion collisions since 
usually only a part of the system - typically around mid-rapidity - is con- 
sidered. Thus energy and conserved quantum numbers may be exchanged 
with the rest of the system, which serves as a heat-bath. There are, however, 
important exceptions when the number of conserved quanta is small. In this 
case an explicit treatment of these co nser ved charges is required, leading to 
a canonical description of the system EH1 and to significant modifications of 
the fluctuations, as we shall discuss below. In the following, we first discuss 
fluctuations based on a grand canonical ensemble and then later point out 
the differences to a canonical treatment. 



2.1. Fluctuations in a grand canonical ensemble 



Assuming we are dealing with a system with one conserved quantum num- 
ber (such as the electric charge, baryon number etc.) the grand canonical 
partition function is given by a 



Z 



E 



Tr[exp(-/3(7J-/iQ))] (1) 



eM-P(H-nQ)) 

states i 

where /3 = 1/T represents the inverse temperature of the system, Q is the 
conserved charge under consideration and fi is the corresponding chemical 
potential. Here the sum covers a complete set of (many particle) states. 
The relevant free energy F is related to the partition function via 

F = -T\ogZ (2) 

We can further introduce the statistical operator 
1 

exp 



PG 



■ exp 



-j3(H - uQ) 



-j3{H-uQ-F) 



(3) 



a Here we restrict ourselves to one conserved charge. Of course the are several conserved 
quantum numbers for a heavy ion collision. The extension of the formalism to multiple 
conserved charges is straightforward. 
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and the moments of the grand-canonical distribution: 

(X n ) G =Tr[X«p G ] (4) 

For a thermal system, typical fluctuations are Gaussian^l anc j are char- 
acterized by the variance defined by 

(SX 2 ) EE (X 2 ) - (Xf (5) 

In the case of grand-canonical ensemble, fluctuations of quantities which 
characterize the thermal system, such as the energy or the conserved 
charges, can be expressed in terms of appropriate derivatives of the parti- 
tion function. Of special interest in the context of heavy ion collisions are 
energy /temperature fluctuations, which are often related to the fluctuations 
of the transverse momentum as well as electric charge/baryon number fluc- 
tuations. 



2.1.1. Fluctuations of the energy and of the conserved charges 

As already pointed out at the beginning of this section when analyzing the 
system created in a heavy ion collisions, one usually studies only a small 
subsystem around mid-rapidity. In a statistical framework, this situation 
is best represented by a grand canonical ensemble, where the exchange 
of conserved quantum number with the rest of the system is taken into 
account. The equilibrium state is then characterized by the appropriate 
conjugate variables, namely the temperature and the chemical potentials 
for the energy and the conserved quantities, respectively. 

As a consequence, energy as well as the conserved charges may fluctuate 
and the size of fluctuations reveals additional properties of the matter, the 
so called susceptibilities, which also characterize the response of the system 
to external forces. 

For example, the fluctuation of the conserved charge in the subsystem 
under consideration is given by 

Similarly, fluctuations of the energy can be expressed in terms of derivatives 
of the partition function with respect to the temperature 

(^ 2 ) G = ^ g Z = -T^F = T^C v (7) 
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Note that the energy fluctuations are proportional to the heat capacity of 
the system, Cy. Thus one would expect that these fluctuations obtain a 
maximum as the system moves through the QCD phase-transition, which, 
among others, is characterized by a maximum of the heat capacity. 

An alternative way ISZlSI i§ to study temperature fluctuations in micro- 
canonical ensemble which are inversely proportional to the heat capacity 

m 

<«">mc = %■ (8) 

However, for situation at hand, which is best described by a grand canonical 
ensemble, energy fluctuations appear to be the more appropriate observable. 

The second derivatives of the free energy, which characterize the fluctu- 
ations, are usually referred to as susceptibilities. Thus we have the charge 
susceptibility 

1 d 2 



and the "energy susceptibility" 



which is related to the specific heat. And just as the specific heat determines 
the response of the (sub)system to a change of temperature the charge 
susceptibility characterizes the response to a change of chemical potential. 
In case of electric charge, this would be the response to an external electric 
field. 

In thermal field theory these susceptibilities are given by correlation 
functions of the appropriate operators. For example the charge susceptibil- 
ity is given by the space-like limit of the st atic time- tim e component of the 
electromagnetic current-current correlator | 41 | 43 | 48 | 22 | 

XQ = -n o(^ = o, g ^o) (ii) 

with 

n„„(fc) = i J d 4 xcM-ikx) (pT*(j u (x)j u (0))) (12) 

It is interesting to note, t hat the charge susceptibility is directly pro- 
portional to the electric mass 1^11 

m l\ = e2 XQ (13) 
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Equation Eq. lllll allows to calculate the electric mass i n any given model 
(see e.g.l 23 l 22 l 48 l 14 l 13 l) and in particular in Lattice QCdEHESI. Since dilep- 
ton and photon production rates are given in terms of the imaginary part 
of the same current-current correlation function - taken at different values 
of u> and q - model calculations for these processes will also give predictions 
for the charge susceptibility, which then can be compared with lattice QCD 
results. As shown in ESI an extraction of \Q from Dilepton data via disper- 
sion relations, however is not possible. For that one needs also information 
for the space-like part of I±oo which is not easily accessible by heavy ion 
experiments. 

Charge fluctuations are of particular interest to heavy ion collisions, 
since they provide a sign ature for the existence of a de-confined Quark 
Gluon Plasma phase I v I vO I Let us, therefore, discuss charge fluctuations in 
more detail. 

Consider a classical ideal gas of positively and negatively charged par- 
ticles of charge ±q. The fluctuations of the total charge contained in a 
subsystem of N particles is then given by 

(SQ 2 ) = q 2 ((5N + - 6N-) 2 ) 

= q 2 [(SNl) + (5N 2 _)-2(5N-SN + )]. (14) 

Since correlations are absent in an ideal gas, (8N_5N + ) = 0. Furthermore, 
for a classical ideal gas, 

(W 2 )„ B = Wmb (15) 
where the subscript MB stands for Maxwcll-Boltzmann. This implies 

(5Q 2 ) MB = q 2 (N+ + AMmb = I 2 <^ch) MB (16) 

where iV c h = N + +N- denotes the total number of charged particles. Taking 
quantum statistics into account modifies the results som ewhat, since the 
number fluctuations are not Poisson anymore (see e.g. I^j) 

<<«V 2 ) G = (N) G (1 ± (n±) G ) = w± (N) G . (17) 

where 

_ fd 3 pn±(p) 2 
(n±> G = fd* P n ± (p) (18) 

Here, (+) refers to Bosons and (— ) to Fermions, and n±(p) represents the 
respective single particle distribution functions. For the temperatures and 
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densities reached in heavy ion collisions, however, the corrections due to 
quantum statistics are small. For a pion gas at temperature T = 170 MeV 
^ = 1.13m 

Obviously, charge fluctuations are sensitive to the square of the charges 
of the particles in the gas. This can be utilized to distinguish a Quark 
Gluon Plasma, which contains particles of fractional charge, from a hadron 
gas where the particles carry unit charge. Charge fluctuations per parti- 
cle should be smaller in a Quark Gluon Plasma than in a hadron gas. 
The appropriate observable to study is then the charge fluctuations per 
entropy. To illustrate this point, let us consider a non-interacting pion gas 
and Quark-Gluon gas in the classical approximation. Corrections due to 
quantum statistics and due to the presence of resonances are discussed in 
detail m EE2ESI j n a neutral pion gas the charge fluctuations are due to 
the charged pions, which are equally abundant 

WmB = Wr+>MB + (^r-)MB (19) 

whereas in a Quark Gluon Plasma the quarks and anti-quarks are respon- 
sible for the charge fluctuations 

«)mb = Ql Wmb + Qd Wmb = I Wmb (20) 

where N q = N u = Nd denotes the number of quarks and anti-quarks. For 
a classical ideal gas of massless particles the entropy is given by 

S MB = 4(iV) MB (21) 

For a pion gas, this yields 

S„ = 4 ({N X+ ) UB + (AV) M b + (N x o) MB ) (22) 

and for a Quark Gluon Plasma 

Sqgp = 4 «AU mb + (N d ) MB + (N g ) MB ) 

= 4 (2 W mb + Wmb) (23) 

where N g denotes the number of gluons. Therefore, the ratio of charge 
fluctuation per entropy in a pion gas is 



< s «i} M B _ i 



(24) 
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whereas for a 2-flavor Quark Gluon Plasma it is 



Wmb _ 1 
Sqgp 24 



(25) 



Consequently, the charge fluctuations per degree of freedom in a Quark 
Gluon Plasma are a factor of four smaller than that in a pion gas. Hadronic 
resonances, which constitute a considerable fractio n of a hadron gas reduce 
the result for the pion gas by about 30 % | 35 | 34 | i eavm g s till a factor 3 
signal for the existence of the Quark Gluon Plasma. 

The above ratio, (<5Q 2 ) / S can also be calculated using lattice QCdE^I. 
Above the critical temperature, the value obtained from lattice agree rather 
well with that obtained in our simple Quark Gluon Pl asm a model here EH. 
More recent calculations for the charge susceptibility « U| give a somewhat 
smaller value, which would make the observable even more suitable. 

Unfortunately, at present lattice results are not available for this ratio 
below the critical temperature. He re one has to resort to hadronic model 



calculations. This has been done in E2ESl using either a virial expansion, a 
chiral low energy expansion or an explicit diagrammatic calculation. In all 
cases, the ratio is slightly increased due to interactions, thus enhancing the 
signal for the Quark Gluon Plasma. 



2.2. Fluctuations in a canonical ensemble 

As pointed out in the beginning of this section, once the number of con- 
served quanta is small, i.e. of the order of one per event, the grand canonical 
treatment, where charges are conserved only on the average, is not adequate 
anymore. Instead the description needs to ensure that the quantum number 
is conserved explicitly in each event. Since the deposited energy is still large 
and is distributed over many degrees of freedom, the canonical ensemble is 
the ensemble of choice. 

Obviously the fluctuations of the energy is identical to the grand canoni- 
cal ensemble, but fluctuations of particles which carry the conserved charge 
are affected. As an example, let us consider the fluctuations of Kaons in low 
energy heavy ion collisions. At 1-2 AGeV bombarding energy, only very few 
kaons are being produced (Nr) ~ 0. 1 EH1 which makes an explicit treat- 
ment of strangeness conservation necessary. For simplicity, let us consider 
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Z l = W I 7^3 ex P (-/3(i?A - Mb)) = <^> MB (28) 



Kaons and Lambdas/Sigmas b as the only particles carrying strangeness. 
In the canonical ensemble, whe re strangeness is conserved explicitly, the 
partition function is given bjl^l 

oo 
n=0 

where Z® cst 1S the standard (grand canonical) partition for all the other non 
strange particles and Z\ , Z\ are the single particle partition functions for 
Kaons and Lambdas, respectively. In the classical limit, these are given by 

Z ° K = 2V Jwf e M~PE K ) = «> MB (27) 

Here, the degeneracy factor of d — 2, 4 for Kaons and Lambdas, respectively, 
take into account the presence of the K° and E particles. Note that Z° is 
simply the number of particles in the grand canonical ensemble in the limit 
of vanishing strange chemical potential. Given th e above pa rtition function, 
the probability P n to find n "Kaons" is given 

by EMU 

Pn = (29) 

where Iq is the modified Bessel function and 

e = Z K Z° A = (N K ) MB (Nl) MB (30) 

Given the probabilities Eq. I|29[) . one can easily calculate the canonical 
ensemble average 

(N 2 K ) C = e (31) 



so that the second factorial moment 

(N(N-1)) C _ 1 e _1 (N K/C 



b In the following we will denote both Lambdas and Sigmas as Lambdas, but include the 
appropriate degeneracy factor to take the sigmas into account 
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This is to be contrasted with the grand canonical result, which follows in 
the limit of e ^ 1. In this case, 

(N*) G -(N) G = (N) 2 G (33) 

so that 

F 2 G = 1 (34) 

Thus for (N) <C 1 the effect of explicit strangeness conservation reduces 
the second factorial moment by almost a factor of two. This suppression of 
the factorial moment due to explicit charge conservation can be utilized to 
measure the degree of equilibration reached in these collisions (see section 
0). 



2.3. Phase Transitions and Fluctuations 

As already mentioned in the introduction, the QCD phase diagram is ex- 
pected to be rich in structure. Besides the well known and much studied 
transition at zero chemical potential, which is most likely a cross over tran- 
sition, a true first order transition is expected at finite quark number chem- 
ical potential. It has been argued ^Hl that the phase separation line ranges 
from zero temperature and large chemical potential to finite temperature 
and smaller chemical potential where it ends in a critical end-point (E) (see 
FigflJ. Here the transition is of second order. There is also a first attempt 
to determine the position of the critical point in Lattice QCD EH. 

As discussed in EH, the associated massless mode carries the quantum 
numbers of er-meson, i.e. scalar iso-scalar, whereas the pions remain massive 
due to explicit chiral symmetry breaking as a result of finite current quark 
masses. 

Fluctuations are a well known phenomenon in the context of phase tran- 
sitions. In particular, second order phase transitions are accompanied by 
fluctuations of the order param eter at all length scales, leading to phenom- 
ena such as critical opalescence^! However, since the system generated in 
a heavy ion collision expands rather rapidly, critical slowing down, another 
phenomenon associated with a second order phase transition, will prevent 
the long wavelength modes to fully develop. In RefP^these competing effects 
have been estimated and authors arrive at a maximum correlation length 
of about £ ~ 3 fm if the system passes through the critical (end) point of 
the QCD phase diagram. 
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In Ref. the authors argue that if the system freezes out close to the 
critical end-point, the long range correlations introduced by the massless 
er-modes lead to large fluctuations in the pion number at small transverse 
momenta. In the thermodynamic limit, this fluctuations would diverge, but 
in a realistic scenario, where the long wavelength modes do not have time to 
fully develop, the fluctuations are limited by the correlation length. In Rcf.^ 
it is estimated that a correlation length of £ ~ 3 fm will result in ~ 5 — 10% 
increase in fluctuations of the mean transverse momentum, which should be 
observable with present day large acceptance detectors such as STAR and 
NA49. Since the precise position of the critical point is not well known, what 
is needed is a measurement of the excitation function of these fluctuations. 

If the system undergoes a first order phase transition, bubble formation 
may occur. Since each bubble is expected to decay in many particle s this 
leads to large multiplicity fluctuations in a given rapidity interval IHE2I 
Fluctuations of particle ratios, on the other hand, should be reduced due 
to the correlations induced by bubble formation"-*! 
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3. Other fluctuations 

In heavy ion collisions, there are many sources of non-statistical fluctua- 
tions. To extract physically interesting information from the observed fluc- 
tuations, it is crucial that we know the sizes of these other fluctuations. 
In this section, we discuss the geometrical volume fluctuations and fluctua- 
tions coming purely from the initial reactions without further rescatterings 
within the context of wounded nucleon model 21 . 



3.1. Volume fluctuations 

In heavy ion collisions, we have no real control over the impact parameter 
of each collision. This implies that geometric fluctuation is unavoidable in 
the fluctuations of any extensive quantities such as the particle multiplicity 
and energy. Furthermore, since the system created by a heavy ion collision 
is finit e, o ne must also consider the thermal fluctuation of the reaction 
volume^. 

To be more specific, consider the multiplicity N. Writing 

N = pV (35) 

where p is the density and V is the volume, the fluctuation of N can be 
written as 

<^ 2 )ebe = <*Abe Webe + (p)ebe (^)ebe (36) 

where the subscript 'ebe' indicates the event-by-event average measured in 
an experiment. Normally, what we are interested in are the fluctuations in 
the thermodynamic limit. In other words, we are only interested in the fluc- 
tuation of the density, (Sp 2 )- The second term containing (SV 2 ) is therefore 
undesirable. 

To make a rough estimate of the volume fluctuation effect, we first 
decompose 

<^ 2 >cbc - <^ 2 ) tll + (^ 2 ) goom (37) 
Assume the ideal gas law PV = NT and using 

<«">. --*•(&), m 

(see 01) the thermodynamic volume fluctuation can be estimated as 

(6V\ h = (Vf/(N) th = (V)/( P ) th (39) 
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or 

(p)l(W)l = (N) th (40) 

which is at least as big as the fluctuation due to the thermal density fluc- 
tuation (c.f. Eq.jnj). 

On top of this, we have purely geometrical fluctuation due to the impact 
parameter variation. To have an estimate, let us simplify the nucleus as a 
cylinder with a radius R and the length L. Assuming geometrically random 
choice of impact parameters, one can show that 

(SV 2 ) oh 2 IRh 3 

\ /gcom _ jOmax ID "max , / ni\ ( A ^ 

~1V) 2 ~ 9^ + + ° {K ™ /R ] (41) 

\ /geom 

if the impact parameter varies between 6 = and b = £> max . The expansion 
is of course valid when 6 max <C R. 

As an example, consider 6 % most central collision of two gold ions. For 
a gold ion, R sa 7 fm and 6 % corresponds to 6 max ~ 3.0 fin. Plugging in 
these values to the above formula yields, 

<^ 2 )geom«(^com/l70 (42) 

or 

(P) 2 (^ 2 )geom = ( Ar ) 2 /170 (43) 

For (N) > 170 then, the observed multiplicity fluctuation can be over- 
whelmingly due to the impact parameter variation. 

This difficulty can be circumvented if the average of the fluctuating 
quantity, say the electric charge Q, is zero or close to it. This can be readily 
seen if one writes 

5Q = 5q(V) chc + (q) chc SV (44) 
ignoring terms quadratic in 8q and SV. Squaring and averaging yield 

Webe - (^ 2 ) cbc (V) 2 + <*Ocbo (<?>L (45) 

assuming that the fluctuations in the charge density q and the volume V 
are independent. Usually, (<5q 2 ) ~ (N) / (V) 2 and from the above impact 
parameter consideration, (SV 2 ) = (V) 2 jC where C = 0(100). 
Therefore, the second term is negligible if 

(^ 2 )ebe (?>L = (Q)cbc /C « Webe (46) 
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or 

(Q)L/Webe«C (47) 

The STAR detector at RHIC counts about 600 charged particles per unit 
rapidity in central collisions and (Q) c ^ c / (A) cbc ~ 0.1. In that case, 

(Q>ebc/Webc«(0-l) 2 x 600 = 6 (48) 

which is much smaller than C — 0(100). If one is to talk about gross 
features, such error would not matter much. However, if a precision mea- 
surement is required, this has to be taken into account. 

3.2. Fluctuation from initial collisions 

One of the important issues in heavy ion collisions is that of thermaliza- 
tion. Seemingly 'thermal' behaviors of particle spectra are very common 
in particle/nucleus collision experiments including e + e _ collisions and pp, 
pp collisions where one would not expect thermalization among secondary 
particles to occur. This seemingly 'thermal' feature results from the way 
an elementary collision distributes the collision energy among its secondary 
particles. Mathematically speaking, this is analogous to passing from the 
micro-canonical ensemble of free particles to the canonical ensemble of free 
particles. By doing so, one incurs an error. However as the size of the system 
grows the error becomes negligible. 

A way to distinguish simple equi-partition of energy from true ther- 
malization through multiple scatterings is to consider not only the average 
values but also the fluctuations. For instance, as explained in section[2| the 
fluctuations of charged multiplicity in thermal system is Poisson or 

«)mb = W*>mb (49) 

However, we do know that in elementary collisions such as pp or pp, the 
fluctuation is much stronger due to KNO scaling 

(^ch)ebc « (^h)L (50) 

Questions then arise: How would ( < ^A r ^ h ) cbc behave if there is no scat- 
terings among the secondaries? Will it resemble the thermal case or retain 
the features of KNO scaling? If it turned out to resemble thermal case, how 
can we distinguish? 
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To answer these questions, let us invoke the wounded nucleon model as 
explained in RefP. In this model, the final charged particle multiplicity is 
given by 



N p 

N ch = J^m (51) 



where N p is the number of participating nucleons (participants) in the given 
event and rii is the number of charged particles from each of the partici- 
pants. The assumption here is that the production of charged particles from 
each nucleon is independent and the produced particles do not interact fur- 
ther. In that case, it is not hard to show 

(N ch ) ehe = <i\g eb » N (52) 

and 



Webe _ (SN p )l be ^ , <<5n 2 ) 



(^ch) obo (N p ) chc 

Here (...) N denotes that this average is taken with respect to a single nu- 
cleon. It is independent of the system size. Therefore the second term is 
where KNO makes its appearance. It states that that (<5n 2 ) N = ckno 
with ckno ~ 0.35 for proton-proton collisions (for instance, see RefP21). 
In Ref. , it is argued that since nucleons inside a nucleus are tightly cor- 
related, the probability to have all nucleons in the reaction volume par- 
ticipate is very high. Therefore (£Wp) ebe / (N P ) ehc should be small and 
(^ch) / (-^ch) = 0((n) N ) which grows with the collision energy. 

To compare with a typical heavy ion experiment, however, two more ef- 
fects need to be considered in addition. First, although the fluctuation of the 
number of participants due to the nuclear wavefunction is small, geometrical 
fluctuations due to the variation of impact parameter can introduce more 
substantial fluctuations in N p . Second, a heavy ion experiment usually can 
observe only a portion of the whole phase space. Therefore when calculating 
(n) and (n 2 ) we must fold in a binary distribution with p = (n) A / (n) iu n 
where Ar/ is the detector window. This implies replacing 

(n) N ^p(n)N, (54) 
<fe 2 ) N ^(l-p)p(n) N +p 2 <5n 2 ) N , (55) 

and 
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(5N?) . (6N?) (56) 

V VI cbc \ P ' gcom v ' 

The fluctuation due to the impact parameter variation was estimated in 
the previous section. Putting everything together then yields 

\ ^^ch )cbc I t\t \ ( ^ ^max \ . / -, \ , , \ 

WWNM = -Tj^T = <-^ch) g eom ( 9^2 R 2 J + 0- ~ P) P c KNO (n) N 

(57) 

Here we assumed that tight correlations renders intrinsic participant fluc- 
tuation very small and also used the fact that (-^Vch) geom = P (Aj))„ com ( n )N- 

A few conclusion can be drawn from Ea. (|57J) . First, if p is sufficiently 
small or if the observational window is sufficiently small compared to the 
whole phase space, the fluctuation of the number of charged particle is 
Poisson. However, this has nothing to do with dynamics. Second, if p is 
sufficiently large, then wwnm becomes significantly larger than 1 since 
ckno is about 7 at RHIC energy. Third, the ratio wwnm depends 
linearly on p. The first point is purely statistical in nature. The second 
and third points can be used to test the validity of this model where no 
secondary scattering occurs. 

What would change in this consideration if thermalization through mul- 
tiple scatterings occurs? Suppose that the total number of produced par- 
ticles is still governed by Ea. (|51fl . In other words, only elastic scatterings 
happened to the secondary particles. Further suppose ^5n 2 ) N is still given 
by the KNO scaling result. In this case, the fluctuations of iV c h in the 
full phase space must remain the same as before and Eq. (|53|l remains valid. 
Since binary process does not really care about clusters in this case, Ea. ij5T|) 
also stays valid. Therefore elastic scatterings alone cannot make any differ- 
ence in multiplicity fluctuations. 

Now let us relax the condition and allow inelastic collisions among the 
secondaries. The essential role of inelastic collisions is to convert energy 
into multiplicity and vice versa. Consider a set of events with the same 
number of participants and same amount of energy deposit. According to 
Eq.(JHHJl, the number of initially produced particles (we will call them 'initial 
particles') have a distribution that is much wider than the Poisson distri- 
bution if (n) N 3> 1. This wide distribution implies that the distribution of 
energy per initial particle also has a wide distribution. As argued above, 
elastic collisions cannot change this situation. If inelastic collisions are al- 
lowed, then an event with smaller number of initial particles would tend 
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to produce more particles since collisions in this event have more available 
energy. In this way, energy is re-distributed evenly and the the multiplicity 
distribution gets narrower. This is, of course, the process of equi-partition 
which lies at the heart of thcrmalization. Therefore if thcrmalization does 
occur starting from a certain energy or a certain centrality, w./v ch must show 
a corresponding behavior changing from approximately (n) N > 1 to lower 
values close to 1 as the energy goes higher or collisions get more central. 
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4. Fluctuations and Correlations 

Fluctuations measure the width of the two particle densities^land therefore 
provides additional information than just the averages. To illustrate this 
point, let us consider a variable x{p) which depends on momentum of one 
particle but does not depend on the multiplicity in each event. Let us further 
consider a generic observable in a given event 

N N 
i=l z=l 

where N is the multiplicity of the event. Since x(p) does not depend on N, 
S is an extensive quantity. 

The event averaged moments of this quantity can be expressed as 

M N m N m 

(^)cbo = m E E - E * m «0~M<») (59) 

m=l — l ifc = l 

where m labels the different events and M is their total number. N m is the 
multiplicity of the event labeled by m. 

On the other hand, the moments of an extensive quantity x(p) calculated 
from n-particle inclusive distribution p n (p\, ...,p„) are defined as 



/ 



dpi...dPn/>n(Pli -;Pn)[x(pi)] 1 ..\x(p n )] 
M N m N„ 

M 



h E E - E K(H)] fei ...K( in )] fc - (60) 

m=l i n — l 

where the sums over i\...i n include only the terms for which all indices 
i\...i n are different from each other. 

One sees immediately that (|59H and (|60ll are related. 



(S)ebc = J d PPl (p)x(p) (61) 
( 5 ' 2 )ebo = / d Pl d P2P2(Pl,P2)x(p 1 )x(p 2 ) + \ dppi{p) [x(p)] 2 (62) 



( 5 ' 3 )cbc = / d Pl d P2dP3 P 3(Pl,P2,P3)x(p 1 )x(p 2 )x(p 3 ) 

3 / dp 1 dp 2 p2(pi,P2) x(p 1 )[x(p 2 )] 2 
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dppiip) [x( P )} 3 (63) 

Similar formulas can be derived for higher moments of S. By setting x(p) = 
1, one can also find that the integral over p n gives the n-th factorial moments 

d Pl ...dp n p n (pi,...,p n ) = (N(N-l)---(N-n + l)) chc (64) 

We have thus established the relation between inclusive m easu rements 
and event-by-event fluctuations for single particle observables 1^1 such as 
e.g. the total transverse momentum, particle abundances 1221 e tc. 

For observables concerning different species of particles, we define 

N 

S a (x) = ^x(pi)S ai (65) 

i=l 

where the 5-function picks out the right species from N particles in the 
event. The event average (S a ) ebe is just as before with the replacement of 
Pi(p) ^ Paip) an d / dpPa{p) — (^a) c b c - The average number of pairs is 
then given by 

(S a {x)Sp{x)) chc = — ^ X (Pi) X (Pj) S PJ S <*i ( 66 ) 
i=l j=l 

In terms of the correlation function, this can be rewritten as 

(S a (x)Sp(x)) chc = J dp a dp p p al 3(p a ,pf3)x(p a )x(p f3 ) (67) 
Again by setting x{p) = 1, one also obtains 

dp a dpp p a p{p a ,pp) = (7V Q iV / 3) cbc (68) 

Similar arguments can be constructed for variables which depend on two 
or more particle momenta, such as e.g. the fluctuations of Hanburry-Brown 
Twiss (HBT) two particle correlations, belong to this class. They are of 
practical interest and will be investigated in heavy ion experiments. Here 
we will restrict the argument to two particles but it ca n b e readily extended 
to multiparticle correlations. The details are given in El 

To summarize, E-by-E fluctuations of any (multiparticle) observable can 
be re-expressed in terms of inclusive multiparticle distribution. In case of 
Gaussian fluctuations the multiparticle distributions need to be known up 
to twice the order of the observable under consideration. 
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Following these straightforward considerations, it is natural to discuss 
fluctuations in terms of two parti cle den sities, or equivalently in terms of 
two particle correlations functions EZEH 

C a p(Pl,P2) = Pap(Pl,P2) - P a {P2)P(j{Pl) (69) 

Here the labels a and (3 denote general particle properties including different 
quantum numbers. Thus correlations between different particle types can 
also be discussed in the same framework. 

Instead of the correlation function C often the reduced correlation func- 
tion 

/_ x _ C a p( Pl ,p 2 ) _ p a p(pi,P2) , . 

n a /3(pi,P2) = — -, — ; — -. — r- — t — ; — -, — r-1 (7U) 

Pa{Pl)Pf3{P2) P a {Pl)Pf3{P2j 

is used. The advantage of using R a p over C a p is that the trivial scaling with 
the square of the number of particles is removed and, therefore, the true 
strength of t he co rrelations is exposed more transparently. In particular, as 
shown in Ref.ESl, the observable R a p has the advantage that it is robust, 
i.e. it is independent of the detector efficiencies to leading order. Also, the 
analysis of elementary reactions such as proton-proto n h as traditionally 
been done based on the reduced correlation function i?o,^!^(see also section 

ED - 

In some cases, such as the charge fluctuations (c.f. 12.1.10 . the fluctua- 
tions of the number of particles in a given region of momentum space is 
considered. In this case one deals with the integrated quantities 

( jv «)a»= / dPocPiipo) (71) 

JAt] 



(N a N ) A - 5 a p (N a ) A = / dp a dp p p 2 (p a ,P{}) (72) 

J An 

where the notation (. . is always to be understood as the event-by-event 
average in a given momentum space region Ary 

1 N m 

i—l 

where 9(p € Arj) — 1 if p falls inside A?y and zero other wise . The correlations 
are then expressed in terms of the 'robust covariances ^1 

p (N a N ) Ar] - (N a ) Av S a0 - (N a ) Av (N p ) Ari 
aP ~ (N a ) Afl (Np) Ar) 
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f Av dpidp2Pap(pi,P2) 



- 1 



(74) 



Sa v d P\p a {pi) J Ar , dp 2 P(3{p2) 



which shares the same virtue as R a f3- 

4.1. Correlations and charge fluctuations 

We will now concentrate on charge fluctuations to illustrate the above gen- 
eral considerations. Here we will mainly work out the relationship between 
the charge fluctuations and the charged particle correlation functions. In 
full phase space, a conserved charge does not fluctuate, or (<5Q 2 ) fuU = 0. 
What we are interested in is, however, fluctuations of the net charge in a 
small phase space window where the effect of this overall charge conserva- 
tion is small. At the same time, this window should be big enough in order 
not to lose information on the widths of the correlation functions. 

From previous section 12. 1.11 we know that if a QGP forms, the charge 
fluctuation per entropy becomes a factor of 2 to 3 smaller than that of the 
hadronic gas. What we are interested in this section is how that is related 
to the properties of the charge correlation functions. Since the net charge 
is Q = N+ — iV_ , the variance is 



Written this way, it is clear that to have a small charge fluctuation, we must 
have a positive correlation between the unlike-sign pairs and it must have 
enough strength to compensate the independent fluctuations of N±. The 
purpose of this section is to elaborate on these points. 

To write the charge fluctuations using the correlation functions, we first 
decompose 



(SQ 2 ) 



(6Nl) An + (8Nl) 



2(5N + 5N-) 



(75) 




(76) 



- 2 



+ / dp a dppC ++ (p a ,p^) + / dpadppC — {p a ,pp) 
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(77) 

The first two terms in the right hand side comes from the fact that integrat- 
ing over like-particle correlations give (N(N — 1)) A „. If {N+) A ~ , 
this can be also rewritten in terms of the robust covariances 

Wa, „, 1 , Wa, v (78] 

V^~ 1 + —^ V ** (?8) 

where 

z/ dyn = R ++ + - 2i?+_ (79) 

and = N + + N- . This is the form advocated in RefP. In this paper, 
the authors argued that in R a p the detector efficiencies cancel out while 
(-^ c h)Ai ? still depends on the efficiency which in modern detectors can in- 
troduce 10 ~ 20% error. 

We know that the single particle distribution p±(p) is proportional to 
the probability density function for single particle momentum. To give the 
two particle distributions similar meaning in terms of underlying correla- 
tions, let us consider a simple toy model. Consider a gas made up of only 
three species of resonances, one neutral, one positively charged and one 
negatively charged. No thermal pions are present. Furthermore let's fur- 
ther assume that the neutral resonance decays into a pair of tt + tt~ and 
when a charged resonance decays it emits only one charged pions. Since 
our model does not contain thermal pions, all final state pions arise from 
resonance decay. We also assume that there is no correlation between the 
resonances themselves. The probability to have a given set of final state 
momentum {{p^}, f° r cnar g e d pions is then given by the product of 

the 7r + and ir~ momentum distribution from each resonance decay: 

n{pt}M}\{ p i}i M ^ M +' M -) 

M M + M- 

= LI Fo(pt,Pa\Pa) II F +(pt\ P b) II F ^Pc\Pc) (80) 
a=l 6=1 c=l 

where {Pi} are the momenta of the resonances just before the decay and 
Mo, M_|_, M_ are the number of the neutral and the charged resonances. The 
F's are normalized to 1. Since we are not so much interested in the distri- 
bution of resonance momenta, we integrate over it with a suitable weight 
(for instance thermal weight) before any other calculation. To simplify, we 
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make a further assumption that the resonance momenta are uncorrelated. 
In that case, c 

M M + M- 

v({ P +},{p7}\mo,m + ,m^ = n f (p+, P -) n u(p+) n /_&,-) pi) 

a=l b=l c=l 

where 

/o(p+ pr) - / dPF (p+,pr\P)P (P) (82) 
/±(pf ) = / dPF±(pf\P)P±(P) (83) 

Here Pi{P) is the momentum distribution for the resonance species «. 

The function V({pf},{p~}\M ,M + ,M_) contains full information 
about the momentum distribution of the final state charged pions given 
the number of underlying resonances. From this distribution all other cor- 
relation functions can be calculated. For instance, the single particle distri- 
bution for positively charged particles is 



M + +M 

(P)= E h d P\ E S(p-Pi)V({pt},{p;}\Mo,M + ,M_) 



x V(M ,M+,M^) 
= {M )h+(p) + {M + )U(p) (84) 

where [dp] is a short hand for integration over all p, and we defined 

h±( P ±) = I dp T f ( P+ ,P-) (85) 



Here V(Mo, M+, M_) is the probability to have the number configuration 
(Mo, M + , M_) in the whole event set and (M±p) is the event average of the 
resonance multiplicities in the full momentum space. Ea. (|84(l thus simply 
states that the number of positively charged particles is given by the number 
of resonances which decay in positively charged pions time the probability 
to for these pions to have the decay momentum p. 



c To be quantum mechanically correct, we need to sum over all permutations of 

and {p~ } and divide by (-M+ + Mo)\(M— + Mq)\. However, since this does not change 

the outcome of our consideration, we will omit that here. 
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Likewise, two particle distributions are 

C±±(p 1 ,p 2 ) = (6M 6M±) h±(p 1 )f±(p 2 ) + (6M 5M±}h±(p 2 )f±(p 1 ) 
+ (SMg) f±( Pl )f±(p 2 ) + (5Mg) h±( Pl )h±(p 2 ) 

- (M±) f±(pi)f±( P2 ) - (Mo) h±( Pl )h±(p 2 ) 

(86) 

and 

C+-(pi,p 2 ) = (SM SM_) h + ( Pl )f_(p 2 ) + (SM Q 5M + ) h_(p 2 )f + ( Pl ) 
+ (6M+5M-) /+(pi)/-(p2) + (<5M 2 ) h+( Pl )h.(p 2 ) 

- (Mo) h + ( Pl )h.( P2 ) + (M ) fo(pi,P2) (87) 

To simplify our consideration, let us regard the charged resonances to 
be iso-spin partners so that we have 

u(p) = f-(p) = m (88) 

For the neutral resonances, /o should satisfy fo{p\,p 2 ) = fo{p 2 ,Pi) which 
leads to 

h+{p) = h-(p) = h{p) (89) 
The average multiplicities are then 

(N±) Av = (M+) [ dpf(p) + {M ) [ dph(p) (90) 

J Art iAjj 

and the charge fluctuation is 

(SQ 2 ) Av = (N ch ) An 

+ [({5M+ - <5M_) 2 ) - (M+) - (M_)] ^ dpf{p) 
-2 (M ) / dp a d Pl3 f o {p a ,p ) (91) 

J An 

We can now consider two situations. First, consider the case where the 
underlying system is a thermal gas of free resonances in the grand canon- 
ical ensemble, where charge is only conserved on the average. In this case 
(<>Mf} MB = (Mj) MB and (6M + SM-) MB = (using Boltzmann statistics 
for simplicity) and the above reduces to 

(5Q 2 )TJ m ' = ( JV ch) A7) - 2 (M ) / dpi d P2 f (pi,p 2 ) (92) 
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If (SQ 2 )^™' is to be substantially different from (N c h) Ari , we need to have 

(M ) ~ (M±) and f Av dpi dp 2 fo(Pi,P2) ~ f Av dpf{p)- In other words > the 
number of neutral resonances have to be large and the correlation sharp. 
In the full phase space, Eq. l|92|l leads to 

/ c-^9\ therm. ,,,,,, 

Wfuu = + (M_) 

(^ch)fou rniaJ (M+) + (M_) + 2(M ) 1 J 

which corresponds to the result obtained in Refill for the thermal resonance 
gas. To see how the finite phase space result <|92[1 differs from the result Q93p. 
consider the extreme case of a flat distribution in the pair rapidity and a 
delta function in the relative rapidity 

^ \ = ^max + y)%max ~ 2/) /g^ 

2y 

max 

/o(yi,y 2 ) - %i - y 2 ) /((yi + y 2 )/2) (95) 

Here we specified our momentum space variable to be the rapidity y and 
is the rapidity of the beam particles in the center of mass frame. In 
this extreme case, it is easy to see that 

<^ch) a"'" = (M+ + M_ + 2M ) p (96) 

and 

(^)7 = (M + + M_)p (97) 
where p — Ary/2y max . Therefore in this case of infinitely sharp correlation, 



______ _ 'full 

/at \ therm. / t\ 7" \ therm. 

{ N ch) Ar , (^ch) fu n 



(98) 



for any values of A77. 

With more realistic / and fo, the relation (|98|l does not hold strictly. 
However as long as 

A?7 > 0^ (99) 

where <r yrel is the width of /o in the y ro i = yi — y 2 direction and the single 
particle rapidity distribution is relatively flat within Afy, Ea. (|98|l should 
hold approximately as long as p is not too close to 1. In this way, one can 
say that the charge fluctuations per charged degree of freedom measured 
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in a restricted rapidity window is a good approximation of the full thermal 
result. From proton-proton collision experiments, we can estimate 

a ytxl « 1 (100) 

Therefore, our rapidity window must be at least that big. 

As noted previously the above results are for the grand- canonical ensem- 
ble. In real-life heavy ion collisions, the concept of grand-canonical ensemble 
cannot be applied to the full phase space since overall charge conservation 
strictly requires (<5Q 2 ) = 0. However, there is no dynamical information 
in this fact. In particular it has nothing to say about whether thermal 
equilibrium has been established within the system. 

To say something about the thermal equilibrium, we need to carve out 
a small enough sub-system and then use the concept of grand-canonical 
ensemble on it. It would be ideal if we could define a fluid cell within 
the evolving fireball that resulted from the collision of two heavy ions. 
Unfortunately, this is impossible for we have no information on the positions 
of the particles, produced or otherwise. Best we can do is to have (pseudo- 
rapidity slices which are supposed to be tightly correlated with the spatial 
coordinates in the beam direction. As shown above, if the system is grand- 
canonical, then this is not a big problem. However, since our underlying 
system is clearly not, we have to know what the effect of having a finite 
rapidity window is and what the measurements in that restricted window 
really signify. In our opinion, what one should try to get is the right hand- 
side of Ea. H93JI which we know is the right grand-canonical limit for the 
resonance gas. 

To take the finite size effect and charge conservation effect into account, 
we make an restriction 5M + — 8M- or M + — M_ = Q c is constant. The 
expression (1911 now reduces to 

(6Q 2 ) Ar] = (N ch ) An - «M+) + (Af_» (j^ */&>)) 

-2 (M ) f d P+ d P - f (p+,p-) (101) 

which differs from the thermal result 1)92(1 by the second term. In the full 
phase space, this expression results in zero as it should. That also means 
that we need to find a way to extract the right-hand-side of Eq. ((93JI from 
the above expression. It will be ideal if we can just measure the second 
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term in Ea. i|l(Jl|l and subtract it. However, to do so we need to know (M±) 
which is not readily available. 

In the literature, two ways of compensating the charge conservation 
effect have been proposed. A multiplicative correction method was proposed 
in RefP^ by the pre sent authors and an additive correction method was 
advocated in Ref!^l To see how the multiplicative corrections work, let us 
assume for simplicity that f. / = f A h = p and rewrite Eq. (|101() as 

(<5Q 2 ) A?; = (1 -p) (N ch ) Av - 2 (Mo) Qf d P+ jf dp_ f (p+,p-) - P 2 ^j 



(102) 



The multiplicatively corrected charge fluctuation is given by 
(SQ 2 ) 



, 9 \ mult. 

(SQ 2 ) 



Ar, (1-p) 

(103) 

This is the origin of the so called "(1— p) correction'^! With the extremely 
sharp correlation function l|95|l and the flat dN/dy Q94|). this yields 

(SQ 2 )2T = (N c *} Av -2{Mo)p={M + +M_)p (104) 
same as Ea. (|97|l . In that case, 



/tit v mult. /tit \ therm. 

{^ch) Av (^ch)f uU 



(105) 



holds for any Ar/. 

With a more realistic / and /o, Ea. ll05|l does not hold strictly. However, 
again as long the condition (|99|1 holds and the single particle rapidity dis- 
tribution is relatively flat within Ar?, Ea. (|lt)5(l should hold approximately 
unless p is very close to 1. 

The authors of RefP^ advocated an add itive method. This method is 
best explained using the expression H78JI • RefPS proposes that this expres- 
sion be corrected by replacing Vdyn — * ^dyn + 4/ (-Nch)fuii- The reason behind 
this correction is that —4/ (-/V cn ) full is the value of Vdyn when N + — N- is 
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fixed and there is no correlation in the system. In our language, this cor- 
rection corresponds to 

= (N ch ) Aii - 2 (Mo) Qf dp+ dp- f (p+,p-) - P 2 )(106) 

Again with the flat f(y) jMJ) and the infinitely sharp f Ea. (|106() 
yields 

mt _ wsr (Mo) ^ 



(iVch) A, (iV ch )™' V Ifcnax / V (M+ + M_ + 2M ) 

where we used p = A?7/2y max . The second term is negligible only in the 
p — > limit. With a more realistic / and /o, this limit becomes 

(108) 

which seems to be of more restrictive use for our purpose of extracting the 
right-hand-side of Ea. il§5|) (or the first term in Ea. (11071) 1. In the full phase 
space limit, 

fu " = 1 (109) 

independent of the specific choice for f . Thus this method overcompensates 
the charge conservation effect slightly. 

One case where we do not need such corrections is when (M±) = or 
(M±) <§C (Mo). This is the case when all particles are produ ced in neu- 
tral clusters. Exactly this type of situation was analyzed in Ref.^. In this 
reference, however, "(1 —p)" correctional was made which actually overes- 
timated the charge fluctuation. How then do we know the relative strength 
of (M±) and (M )? 

From Ea. i|101|l . it is clear that the unwanted second term scales like 
p 2 . The first term (N c h) of course scales like p. The last term in Ea. (|101|) 
however scales differently as the size of Ary varies. If A77 is much smaller than 
the correlation length, this term varies like p 2 . However, if Arj exceeds the 
correlation length, then integration along p lc i = p + — p- ceases to change 
and hence it varies like p. This observation suggests the following method: 
Vary the observation window size Ar/ and plot {3Q 2 ) Ari as a function of 
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p = (N c h) A , n I (-Nch)total- According to Ea. (|101|) . this plot can be described 
by a quadratic polynomial in p 

{6Q 2 ) An = ap-bpi (110) 

If the single particle distribution is flat (as is the case for rapidity distribu- 
tion), then p ex Ar] and the parameterization (<5Q 2 ) Aj) = a' Arj — b' (Ar/) 2 
can be used as well. 

The coefficients a and b are not truly constants but vary slowly as a 
function of p. For small p, 

a « (iVch) ful] = (M+ + M_ + 2M ) (111) 

After Ar] exceeds the correlation length in /o, the /o term becomes linear 
in p and the coefficient b becomes. 

b w (M + + M_) (112) 

Therefore by measuring the coefficients in these limits, one can estimate 
the relative strength of (M±) and (Mo) and fully correct Ea. i|lUl[l . 

At this point, one should ask how all this can be reconciled with the 
thermal fluctuations we considered in previous sections. The results are 
certainly similar. However, it seems that we have used no thermal features 
at all in deriving Eas. H101ll03fl . The expression (|102[l rather suggests that 
the these expressions result from binary process where a particle has a 
probability of p to end up in the detector. Of course, when N is large and 
p is small in such a way to have Np fixed, binary process becomes Poisson 
process but that does not mean that the underlying particle energies are 
thermally distributed. So how can the results of this section be reconciled 
with the QGP result from previous section? 

Recall that the arguments in section 12. 1.11 that led to our conclusion of 
reduced charge fluctuation depended on two facts. One, the single particle 
distributions in particular the particle abundances are thermal. Two, the 
number fluctuations are all approximately Poisson. As Ea. 1)84(1 shows, the 
single particle distributions in this section do depend on the underlying 
momentum distribution. In particular, if one is to argue relationships such 
as Ea. 1)21(1 hold, one must have thermal momentum spectra. On the other 
hand, it doesn't really matter to the rest of the arguments whether the 
Poisson nature of the number fluctuation fEas. ((19120(1 is a result of having 
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a thermal system or just a result of having a small detector window. There- 
fore, the fact that Eas- HlUll ll(J2[l are not the consequences of underlying 
thermal system does not negate our previous conclusion. 

Although the argument given above does reconcile the thermal consid- 
eration and the results from present section, it is still unsatisfactory in 
some aspects. The essential point here is our inability to have a sub-system 
in the sense of grand canonical ensemble. This is because we can't carve 
out a sub-system in the coordinate space to observe. Unfortunately this 
limitation is unavoidable. The best we can do is to define sub-systems in 
(pseudo-)rapidity space and argue that strong longitudinal flow strongly 
correlates the (pseudo-)rapidity and the coordinate in the beam direction. 
It would be desirable to work-out ensemble average taking into account this 
fact. That analysis, though, is still to be carried out. 

Another question to ask at this point is whether the underlying corre- 
lation /o(pi,P2) really corresponds to resonances or simply indicates that 
when the pions hadronize from a QGP most of the times they are pair 
produced. One way to distinguish the two scenario may be to actually mea- 
sure the correlation length. Since the unstable neutral meson masses are 
all larger than twice the pion mass, thermal resonance gas must exhibit 
a certain characteristic momentum correlation length between the charged 
pions. For instance, if a p° at rest decays into two pions, they are 3.5 units 
of rapidity apart along the line of the decay. Pions directly coming out of 
hadronizing QGP on the other hand have no such strict constraint imposed 
on their correlation and should exhibit much sharper correlation. A lack of 
good hadronization scheme from a QGP makes this consideration difficult 
to substantiate. However, these points need to be further clarified. 



4.2. Correlations and Balance Function 

A balance function is a particular way of combining the correlation func- 
tions to articulate an aspect of the system. The balance function con- 
structed by Pratt et.al. uses the number of like-sign pairs and unlike-sign 
pairs within a given phase space volume Ar/: 

B( V \A V ) 

= 1 r (JV + _(7 ? |A7?)) (N- + ( V \Ar,)) _ (N ++ ( V \A V )) _ (N—(r,\ At,)) 
2 [ (N-(An)) + {N + (Arj)) {N + (A V )) (N-(A V )) 

(113) 
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For instance, in this expression N-\ (rj\Ar]) is the number of unlike-sign 

pairs which are rj apart from each other within the window A77. To relate 
(Nij) to to correlation functions, first we express the number of pairs within 
Arj and rj apart as 



N^Ar,) = £ 6(p? e A v )6(p" e A V )6(\p? - Pil - v) (114) 

ij 

Integrating over the expression (|81|l give 



(Nab^lAri)) = / dp 1 dp 2 5(\pi - p 2 \-r)) p ab (p 1 ,p 2 ) (115) 
J An 

With our model of resonance gas, the balance function can be expressed as 
B^A-q) » 1 [2 (Mo) / dpi d P2 S(\ Pl - P2 \ - r,) / (pi,P2) 



Ai) L J A 77 

+ «M + ) + (M_)) / ^1 d P2 tfClpa -p 2 | - t?) /(pi)/(pa) 

(116) 

assuming that (Q) Ar) <C {N c b) A . From the expressions Ea. l|lltj|) and 
Eq. (|101jl it is clear that 



<<5Q 2 > 



A// 



y dr]B(r]\Ar]) (117) 



(^ch) Al) 

Balance function is originally devised to detect the change in the unlike- 
sign correlation function. Just as the charge fluctuation, this is possible only 
if Mo ~ M± and fo is a sharply peaked function of rj. Pratt et.al. argued 
that if a QGP forms, the width of the balance function B(rj\Arf) will be 
reduced by a factor of 2 compared to the width of the resonance gas balance 
function. Recent measurement by STAR collaboration^! indicates that go- 
ing from peripheral collisions to the central collisions there is about 20 % 
reduction in the width of the balance function although it is yet not clear 
whether this is indeed the signal of QGP formation. 
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5. Observable fluctuations 

Although the motivation for studying fluctuations is similar to that in solid 
state physics, the observables in heavy ion collisions are restricted to corre- 
lations in momenta and quantum numbers of the observed particles. Spatial 
correlation are only indirectly accessible via Fourier transforms of momen- 
tum space correlatio ns, and thus limited. An example is the Bose Einstein 
correlations (see e.g.'"*'). 

Consequently, the fluctuations accessible in heavy ion collisions are all 
combinations of (many) particle correlation functions in momentum space. 
In addition, as discussed in the previous section, even for tight centrality 
cuts there are fluctuations in the impact parameter which may mask the 
fluctuations of interest. In the thermal language, these impact parameter 
fluctuations correspond to volume fluctuations. Consequently, one should 
study so called intensive variables, i.e. variable which do not scale with the 
volume, such as temperature, energy density etc. 

Another issue is the presence of statistical fluctuations due to the finite 
number of particles observed. These need to be subtracted in order to access 
the dynamical fluctuations of the system. 

Finally, although this is outside the expertise of the authors, there are 
fluctuations induced by the measurement/detector, which also contribute 
to the signal. Those need to be understood and removed/subtracted as well. 

Let us return to the first two issues. As already discussed in some detail 
in section 0] the fluctuations can always be described in terms of the in- 
clusive many-particle densities in momentum space p n (pi, ■ ■ ■ ,Pn)- For the 
simplicity of the discussion let us concentrate on two particle correlations, 
which fully characterize fluctuations of Gaussian nature d . 

Let us start the discussion by defining the number of particles with 
quantum numbers a = a%, . . . , a n in the momentum space interval (p,p + 
dp) in a given event 3 



dN a 

_gent ( n8 ) 



d In case of fluctuations of positive definite quantities, such as transver se m omentum or 
energy,the appropriate distribution is a so called Gamma- distribution EH1 
e In this section, we will use the term 'event' and 'member of the given ensemble' inter- 
changeably. We also use 'event average' and 'ensemble average' interchangeably. 
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and its fluctuations 

Sn a p =n°- <n°) (119) 

In the case of ideal gas, (rip) takes the form of the Bose- Einstein distribu- 
tion or the Fermi-Dirac distribution depending on the spin of the particles. 
The mean value of an observable f 

x = J2 x >p ( 12 °) 

is obtained by averaging over all the events in the ensemble 

(*)=5>£K> (121) 

p . a 

Note that X defined in this way is an extensive observable. 

The goal of studying fluctuations is to see the effect of dynamics in 
terms of non-trivial correlation. If the system under study is totally devoid 
of correlation, then the single distribution function alone must be able to 
describe all the moments of X since higher order correlation functions would 
be just products of single particle distribution functions. Any deviation from 
this behavior is a sign of non-trivial correlation. To make a quantitative 
calculation, let us define the single particle probability density function 

^ incl= jh Pa{p) (122) 

where p a (p) is the single particle inclusive distribution function introduced 
in section 2| Note the factor 1/ (N a ) in the definition of n£ mc which nor- 
malizes this distribution to unity. Note also that event average has already 
performed and the only variable left to average over is the single momentum. 
We define the average taken with this probability density as the 'inclusive 
average'. 

^ nCl = E^ inCl (123) 

p,a 

For a single particle observable X 

(124) 



f Here we use sums over momentum states, as appropriate for a finite box. The conversion 
to continuum states is a usual, ~ * ^ J (2 J) 3 ' wnere ^ ' s ^ ne volume of the system. 
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For the discussion of fluctuations the relevant quantity to consider is 

A«f = (5n%5nP) (125) 

where Sup 1 is defined in Ea. (|119fl . For an ideal gas, only the identical par- 
ticles are correlated and ((Srip) 2 ) = ujp (rip). Therefore, 




(126) 



Here 

w« = 1 (127) 

for a classical ideal gas, and 

W «=(l±(n£» (128) 

for a boson (+) and a fermion (— ) gas. Before we continue, let us step 
back and try to understand eqs. (|126[) to l|128|) in simple terms. The term 
Ili 3ai/3i ensures that only identical particles are correlated. For example, 
7r + and 7r~ , are both pions but with different charge. Therefore they are 
not correlated in an ideal gas, i.e. (Stt + 5it~) = 0. Also there are no effect 
due to quantum statistics, as the particles are different. 

In the presence of dynamical correlations, the basic correlator (|125fl may 
contain off-diagonal elements in momentum space and/or in the space of the 
quantum numbers. For example the presence of a resonance, such as a po, 
correlates the number of tt + and 7r~ in the final state due to resonance decay. 
Also their momenta are correlated. In addition, the occupation number n p 
in a given momentum interval may be changed. This would be for instance 
an effect of hydrodynamical flow. 

Any two- partic le observable can be expressed in terms of the basic cor- 
relator Ap | 12 | 56 | p or g erier j c observable X as given by eq. I|120|) the 
variance is given by 

(sx 2 )= Yl a p:^p< ( 129 ) 

p,q,a,{3 
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To illustrate the formalism, let us consider the net charge 8 

Q = N+ — N- = q a rip (130) 

Here, q a is the charge of the particle whereas all other quantum numbers 
are still being summed over. Likewise the number of charged particle is 
given by 

N ch = N+ + iV_ = k-l n° (131) 

P,CK 

The variances of these quantities are given by 

(SQ 2 ) = (SNl) + (SN 2 ) - 2 (6N+5N-) 

= E( a pY + a p t- 2A p V) ( 132 ) 

p,i 

(5N 2 h ) = (SNl) + ( SN -) + 2 (SN+6N-) 

= E ( A ™ + a pX + 2A ir) ( 133 ) 

Here the basic correlator A+^ + still contains a sum over all other quantum 
numbers, which we have suppressed. 

A+;+= E *#r }A+M (134) 

ai ^charge, fli ^charge 

and similar for the other combinations. This sum simply means that all 
charged particles are included independent of their flavor, spin etc. In the 
following, we will always use this short-hand notation with implicit sum- 
mation over all the quantum numbers not specified. 

Obviously all the information revealed by the fluctuations is encoded 
in the basic correlator and the fluctuations of different observables 

expose different "moments" and elements of the basic correlator. Using 
the notation from the previous section, the basic correlator for a f3 
is nothing but the 2-particle correlation function defined in Eq. I|69|I . For 
the identical particle correlator, Ea. l|125|) differs from Ea. (|69|l by a single 



g Note the difference in notation: whereas 8n p refers to the fluctuations in the momentum 
interval (p,p + dp), 8N refers to the fluctuations of the total (integrated) number of 
particles. 
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particle distribution so that 




p a p(p, q) - p a (p)pp(q) + 8 af3 $pq Pa (p) 
C a p(p, q) + S a/ 3 S pq p a (p) 



(135) 



Therefore, an explicit measurement of the two-particle correlation func- 
tion C a /3 would be extremely valuable, since all (Gaussian) fluctuation in- 
formation can be extracted from it by properly weighted integrals. 

If the system has no genuine two particle correlations, then the basic 
correlator will be 



similar to the ideal gas, only that (n p ) does not have to follow a Boltzmann 
distribution. Using the inclusive single particle spectrum (|122[l as 

momentum distribution, the relation (jlr?6|) may b eut ilized as an estimator 
of the fluctuations due to finite number statistical which follow Poisson 
statistics. A more detailed discussion on this aspect will be given below. 

Any correlations, on the other hand will lead to a non Poisson behav- 
ior of the basic correlator. The ideal quantum gases discussed above are 
examples where the symmetry of the wavefunction introduces 2-particle 
correlations which either reduce (Fermions) or enhance (Bosons) the fluc- 
tuations. Also global conservation laws such as charge, energy etc. will affect 
the fluctuations. However, if only a small subset of the system is discussed 
these constraints will be negligible 11 . As the fraction of accepted particles 
increases, these constraints may become more relevant, however, and thus 
need to be properly accounted for. This has been discussed in some detail 
in the previous section. 

As mentioned in the beginning of this section, in order to avoid con- 
tributions from volume / impact parameter fluctuations it is desirable to 
study so called intensive quantities, i.e. quantities which do not scale with 
the size of the system. Examples which are currently explored experimen- 
tally are the mean transverse momentum and mean energy fluctuations as 
well as fluctuations of particle ratios. 



h This is the same argument which leads to the canonical and grand-canonical description 
of a thermal system, as briefly discussed in section l2l 



(136) 
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5.1. Fluctuations of Ratios 

Let us begin with a general discussion of fluctuations of ratios, since mean 
energy or mean transverse momentum as well as particle ratios are all ratios 
of observables. Consider the ratio of of two observables A and B 1 

B = E 6 X- ( l3? ) 

To find the average and the fluctuation of the ratio 

Rab = 4 ( 138 ) 

we first write A = (A) + 5 A and B = (B) + SB and expand the numerator 
and the denominator to get 



Rab 



(A) (A) (SA SB\ 

Wi + Wi{jA)-Wi) +oi5) (139) 

where 0(S 2 ) indicates terms that are of quadratic and higher orders in 8 A 
and 5B. Since A and B are extensive observables, the neglected terms are 
at most 0(1/ (N)) which we will neglect from now on. From Ea. (|139l) . it is 
easy to see that 

(Rab) = t4t ( 14 °) 



(B) 



and the variance 

?2 



X AB/ (B) 2 \ (A) 2 (B) (A) (B) j 

Using the basic correlator this can be rewritten as 

For an ideal gas, this simplifies to 

, 2 \ ideal _ (A) 2 / (qg) 2 (6g) 2 qgbg \ 



'Defined this way, both A and B scale with N, the number of particles in the final state. 
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Obviously, deviations from the ideal gas value will give us insight into the 
dynamical correlations of the system. In order to expose these deviations, 
very often one compares to the fluctuations based on the 'inclusive' single 
particle distribution. The latter estimates the contribution of finite number 
fluctuations to the observed signal. In our formalism, this means that one 
evaluates equation (|142|l with the fully uncorrelated basic 'correlator' as 
defined in equ. (|136|l . In addition, as it is common in the literature, one 
replaces the event averages (...) by inclusive averages t-t-t , which simply 
mean multiplying by the appropriate factors of (N). This let us define 

cl 



(SR) 2 

(i^ EE ^> ( '^ + jffiL_ 2 ^p_y (144) 

In absence of any dynamical correlations 

((6Rf) = ^W md (145) 

This observation has lead to the measure of dynamical fluctuations Od ynamic 
1^01 which is given by 

dynamic = ((«?) 2 ) " T^W' (146) 

The fi rst such measure to be proposed has been the so called $ variable 
| 27 | 46 | which m terms of our variables here is given by 



$ = y/(N) ((6R) 2 ) - V (SR) 2 . (147) 

The ratio has also been proposed ^1 

(N) ({6R) 2 \ , s 

F ^ 2jLA ' (148) 
{SR) 2 

which has the advantage of being dimensionless. 

Aside from simply subtracting the expected fluctuations from an uncor- 
related system, the sub-event method has been developed in ^4 

We note, that given the same momentum distribution the inclusive re- 
sult is up to a trivial factor identical to that of an classical ideal gas, i.e 
^dynamic = $ = 0orF=lin this case. However, in reality, the inclusive 
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momentum distribution usually differs from a Boltzmann shape due to ad- 
ditional effects such as hydrodynamic flow. Notice also, that the results for 
Bose and Fermi gases already differ from the inclusive estimator, reflect- 
ing the correlations induced by quantum statistics. For a Bose gas, $ > 0, 
F > 1, whereas for a Fermi gas $ < 0, F < 1. For the observable to be 
discussed below, the corrections to F are a few percent ESI 

Finally, the above formalism allows us to discuss more general correla- 
tions between ratios of observables. Lets introduce two more observables C 
and D 

C = £c£n» D = J2d«n%. (149) 
Then, the most general correlation of ratios can be written as 

"i) '(§:■>- 

{A)(C} ( (6A5C) (SB6D) (6 ASP) (6B6C) \ 
(B) (D) { (A) (C) + (B) (D) (A) (D) (B) (C) J [ > 
Using the basic correlator l|125|) we obtain 

(C) 

(A)(D) (B)(C)) ™ [ ' 

Note, that again this correlation function represents nothing but a specific 
moment of the basic correlator. Obviously, the replacements C = A and 
D = B leads to the ratio fluctuations discussed above. This general cor- 
relation will become useful below, when discussing charge dependent and 
charge independent tran sverse momentum fluctuations, as proposed by the 
STAR collaboration E2. 

After having developed the general formalism for fluctuations of ratios 
it is now straight forward to discuss the actual observables. 

5.2. Fluctuations of the mean energy and mean transverse 
momentum 

The mean energy and transverse momentum are defined as 

E _ E P n p E{p) 




c 
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p t E p n P Pt(p) 

where E and (Pt) denote the Energy and transverse momentum of an event 
with N particles in the final state. Obviously, these observables are ratios 
and the above formalism can be readily applied with the substitutions 

A = E; or A = P t 

B = N (153) 

So that for X being either the transverse momentum of the energy 

After rearranging terms one finds for an ideal gas 

— Y 

(N) 2 V 

where X stands either for the energy E or the transverse momentum P t . 

In addition to the transverse momentum fluctuations for all charged 
particles, one can investigate the p t fluctuations of the negative and positive 
charges independently as well as the cross correlation between them. Let 
us define 

Pt =TT = ±/ \ ( 156 ) 

N + E P n P (p) 

where n^{p) is the momentum distribution of the positively (negatively) 
charged particles. Using Ea. (|142|l . we then have 

Pt(p)pt(q) . 1 Pt(p)+Pt(q) 





2 

(155) 




±,± 



(Pt) 2 (N*) 2 (A f± > (Pt ) 

„±\\ /„ f„\ /±\\ A±,± 



v ,, ^ (PM ~ <Pf » (Ml) - (Pt)) (157) 

(^i) p,q 

where to get the last line we used (pf) = (Pt) I (N±)- For the unlike-sign 
pairs, we use Ea. (|15ip to get 

<* ft+ * ft "> = (AT + )V) ^ {Ptiq) " <Pt+)) " {P7)) A ™ {m 
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Therefore it is natural to define the measure of dynamic fluctuation as 



■incl 



A*\=(N ± )((5ptf)-(6ptf (159) 



-inc 



For unlike-sign pairs, Sp^8p t — so that 

AoJ+_ = y/{N+) (JV_) (SptSpT) (160) 

This then allows us to define the 'charge- independent' (CI) and 'charge- 
dependent' ( CD) combinations as used in the STAR and the PHENIX 
collaborations^! 



(N) A^ I>CD (p t ) = (N+) A* 2 + (iV_) Ac 2 ± 2^{N + ) (N_)Aa 2 + . 

(161) 

Here the (+)-sign in front of last term leads to the 'charge-independent' and 
the (— )-sign to the 'charge dependent' combinations, respectively. Under 
the quite reasonable assumption that 

i^l <£> (g±£) = {Pt) (162) 

(N+) (N-) ~ (N+ + A_) ~ {Pt) 1 ' 



and (8pt ) 2 = (<5p t ) 2 = {Spt) 2 , this can be rewritten as 
(N ch ) Aa 2 ChCB (p t ) = 

P,8 



- 2 (iV ch ) (<5p t ) 2 (163) 

where the assignment of plus and minus signs are as in Eq. (|161(l . 

Notice the the combinations of the basic correlator A in the 'charge- 
dependent' combination is the same as that of the net charge fluctuation 
Ea. (|132|) . Likewise, the 'charge-independent' combination resembles the 
fluctuation of the fluctuations of the number of charged particles, eq. (|133l) . 
As we discussed in detail in section 0] charge conservation suppresses the 
fluctuations of the net-charge in a subsystem. Therefore, one would ex- 
pect, that also the 'charge-dependent' transverse momentum fluctuations 
are suppressed, since they are nothing but a different moment of the same 

cor relator. This is seen in the preliminary data of the STAR collaboration 
51 
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Other correlation, such as resonanc es a nd flow may also affect the trans- 
verse momentum fluctuations. In RefP these have been estimated and 
found to be small, of the order of 1 — 2%. 

5.3. Fluctuations of particle ratios 

The fluctuation of particle ratios are also presently studied. Consider the 
ratio of particles "1" and "2" 

Ru = f (164) 
then with A = JVi and B — N2 the dispersion is given by eq. I|141[l as 

W .jtf(W t W-,M) (165) 

Notice that the last term introduces correlations between the particles 
which may reduce the fluctuations. One example for these correlations are 
hadronic resonances. Consider, for example the ratio of positively and neg- 
atively charged pions R = The the presence of neutral rho- mesons, 
which decay into n + and tt~ reduce the fluctuations. This can be easily 
understood by considering a gas made out of only neutral rho mesons. In- 
dependent, how large the fluctuations of these are, in every event on has 
as many n + as n~ , and hence the fluctuations of their ratio vanish. This 
effect can be utilize d to estimate the number of resonance in the system at 
chemical freeze-out 

Next, consider the situation where one particle species is much more 
abundant then the other, i.e. (N 2 ) ^> (Ni). This is the case of the fluctua- 
tions of the kaon to pion ratio, as investigated by the NA49 collaboration. 
Assuming that the number fluctuations are approximately Poisson, then 
the second term in eq. ()165|) dominates 

(AM 2 " W ~ (N,f 1 j 

so that 

/rn2 \ 1 

i-^Z ~ " (1 + O ((AT 1 )/(iV 2 ))) (167) 

In other words, if one particle species is much more abundant then the 
other, the correlations among the particles have to be very strong to be 
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visible in experiment. To give some numbers, in a standard resonance gas 
1^ the correlations due to resonance lead to 30% corrections in case of the 
n + /ir~ ratio and only to 4% in case of the K/ir ratio 

Thus, ratios of equally abundant particles are best suited to expose 
possible correlations. The simplest, and arguably most interesting is that 
of positively over negatively charged particles 

R+- = j± (168) 

\ + 1 \ ( N+ ) 2 ( N _f (N-) (M) J 

In the limit that the net charge (Q) = (N+ — N-) is much smaller than the 
number of charged particles (iV c h) = (N + + N-), or (Q) <C (N c h) 

(B+-)*l; (N + )~{N.)~^- (no) 

so that 

(«*+-> = TTT^ ( SN l + SN - - 26N + 6N ~) 

(171, 

Since the number of charged particles is a measure for the entropy of the 
system, (iV c h) oc S, the observable 

, „ . (SO 2 ) (SO 2 ) 

is a measure for the charge fluctuations per entropy. And this, as discussed 
in detail in section 12.1.11 is an observable for the existence of the Quark 
Gluon Plasma. As discussed in detail in ESI for a pion gas, £> P ion-gas = 4 
whereas for a QGP, -Dqgp — 1 — 1-5, where the uncertainty arises from 
relating the entropy S with the number of charged particles (iV c h) . Hadronic 
resonances introduce additional correlations, which reduce the value of the 
pion gas to -Dhadron-gas — 3, but still a factor of 2 larger then the value for 
the QGP. 

(so 2 \ 

In principle, one could directly measure the ratio W^T > without go- 
ing through ratio fluctuations. However, since the net charge is an exten- 
sive quantity, this will introduce volume fluctuations into the measurement. 
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Only in the limit that the total charge of the system is zero, volume fluc- 
tuations do not contribute to the order considered here . 

The key question of course is, how can these reduced fluctuations be 
observed in the final state which consists of hadrons. Should one not expect 
that the fluctuations will be those of the hadron gas? The reason, why it 
should be possible to see the charge fluctuations of the initial QGP has to 
do with the fact that charge is a conserved quantity. Imagine one measures 
in each event the net charge in a given rapidity interval Ay such that 

Ay coU < Ay < Ay 

max 

(173) 

where Ay max is the width of the total charge distribution and Ay co u is the 
typical rapidity shift due to hadronization and re-scattering in the hadronic 
phase. If, as it is expected, strong longitudinal flow develops already in the 
QGP-phase, the number of charges inside the rapidity window Ay for a 
given event is essentially frozen in. And if Ay 3> Ay co u neither hadroniza- 
tion nor the subsequent collisions in the hadronic phase will be very effective 
to transport charges in and out of this rapidity window. Thus, the E-by-E 
charge-fluctuations measured at the end reflect those of the initial state, 
when the longitudinal flow is developed. Ref. ^1 arrives at the same conclu- 
sion on the basis of a Fokker-Planck type equation describing the relaxation 
of the charge fluctuation in a thermal environment. 

In Fig|21 we show the results of an URQMD calculation ^1 (left fig- 
ure) , where the variable D is plotted versus the size of the rapidity window 
Ay. For large Ay the results have to be corrected for charge conservation 
effects; if all charges are accepted, global charge conservation leads to van- 
ishing fluctuations (open symbols in Fig[21). This can be corrected for as 
explained in the previous section (for details see^j). The resulting values 
for D are shown as full symbols in Fig|2] They agree nicely with the predic- 
tion for the resonance gas, as they should, since the URQMD model does 
not contain any partonic degrees of freedom. For small Ay the correlations 
imposed by the resonances are lost, because only one of the decay products 
is accepted. As a result we see an increase of D. The right figure shows 
a comparison of several transport models, including the Parton Cascade 
model This model starts with partons in the initial state, but has some 
model for hadronization included. If the general ideas about the reduced 
charge fluctuations are correct, this model should lead to smaller values of 
D, and it does. 
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Fig. 2. Left: Charge fluctuations for different rapidity windows as obtained in the 
URQMD model. Open symbols: without correction for global charge conservation. Full 
symbols: With correction for global charge conservation. Right: Results for different 
transport model including the parton cascade 1^21 These results are corrected for charge 
conservation. 



For very small Ay, when (N) ~ 1, the ratio D is not well defined for 
events with iV_ = 0, and, therefore, cannot serve as a observable. Alter- 
native observables, measuring the same quantity have been proposed and 
studied inEHl 



5.4. Kaon Fluctuations 

As discussed in section at low energy the strangeness conservation has 
to be treated explicitly. This can be utilized to determine the degree of 
equilibration reached in these collisions 1251. By measuring the ratio 

F 2 = W"* a -*)> (174) 
(N K f 

of the number of kaon pairs over the square of the inclusive number of 
kaons a factor two sensitivity on the degree of equilibrium can be obtained. 
From transport models (see e.g.^J) it is known, that most of the kaons are 
produced in secondary collisions, i.e. during the evolution of the fireball. 
Transport calculations also find that the equilibration time of kaons is of 
the order of 500fm/c, which is about ten time the lifetime of the system. 
Hence these models would not predict equilibration of the kaons. However, 
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observed particle ratios including kaons are also consistent with a thermal 
model description (for a review see^J). The measurement of F% can resolve 
this issue and point to new physics such as multi-particle processes, in 
medium effects etc., if indeed it is found to be consistent with equilibrium. 
This is demonstrated in Fig|3| where the time evolution of F2 for several 
initial kaon numbers is shown. In all cases, F2 quickly rises close to F% — 1 
before it settles at the final equilibrium value of F^' ~ 1/2. Thus, by 




Fig. 3. Time evolution of F2 for various initial kaon numbers. The thin lines are ana- 
lytical results for early times. 



measuring F% one can directly determine how close to chemical equilibrium 
the system has developed, before it freezes out. If the predictions of the 
transport models are correct, then a value of F2 ~ 1 should be found. If, 
on the other hand, equilibrium is indeed reached in these collisions, then 
F2 — |. In principle a similar measurement can also be done at higher 
energies for charmed mesons. To which extent this is technically feasible is 
another question. 

Let us conclude the section by mentioning other observables. In the 
context of so called DCC production the fluctuations of the fraction 
of neutral pions is considered a useful signal. Also the fluctuations of the 
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elliptic flow has be proposed as an useful observable ^4 which may reveal 
new configurations such as sphalerons, which could be created in heavy ion 
collisions. 
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6. Experimental Situation 

6.1. Fluctuation in elementary collisions 

By now, it is clear that Quantum Chromodynamics is the right theory of 
strong interaction. However, essential non-perturbative nature of the strong 
interaction still remains a difficulty. In heavy ion collisions, this is more 
manifest in the sense that even the hard part of the spectrum is strongly 
influenced by the surrounding soft medium through multiple scatterings. 
Therefore, before we begin to consider the experimental results from heavy 
ion collisions, it is crucial that we understand elementary collision results 
such as proton-proton collision results since these elementary collision re- 
sults can provide unbiased baseline. 

A large amount of data for pp collisions at various beam ener gies were 
taken from Fermi Lab in the i970'sEE3SHEni ( S ee also Ref.EU). Among 
them, data at the be am energy of 205 GeV were most thoroughly analyzed 
by Kafka et.alPSESl First, let us consider the 'charge fluctuation study'. 
In the case of pp collisions, the definition of 'charge fluctuation', however, 
differs from ours. Let us define the charge transfer 

u{y) = \[Qf{y)-Q b {y)\ (175) 

where Q / (y) is the sum of the charges of the particles with rapidity larger 
than y and Qb(y) is the sum of the charges of the particles with rapidity 
smaller than than y. We then define the charge transfer fluctuation 

(^(y) 2 ) cbc = (u 2 (y)) chc -Hy)) 2 chc m 

This quantity is usually referred to as 'charge fluctuation' in literature 
dealing with proton-proton collisions. This charge transfer fluctuation is 
certainly not the same as what we have been discussing so far which is the 
charge fluctuation within a given phase space window Arj. However, as we 
will show shortly, (<5u(z/) 2 ) is intimately related to the charge correlation 
functions. By considering the charge transfer fluctuation, we can then put 
a constraint on possible forms of the correlation functions. 

To do so, consider again our simple model defined by Eq.JSTJ. As before, 
we impose the conditions M + — M_ = Q c = 0, /+ = /_ = / and h+ = 
h- = h. After straightforward but tedious algebra, we obtain 

poo py 

(5u(yf) chc = «M + ) + (M_» / dy' f(y') / dy" f(y") 

J y J —oo 
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r°° rv 

+ 2 (M ) / dy+ / dy-f (y+,y-) (177) 

using the rap idity y as the phase space variable. 
In Ref!"*", it is shown that the data satisfies 

(Su(yf) chc « 0.62^ (178) 

at the beam energy of 205 GeV and the shape of dN c h/dy is very well 
represented by a Gaussian with <jy ss log( v / s/2mAr) w 2.2. This result 
puts a condition on possible forms of / and Jo- In our model, the rapidity 
distribution of the charged particle is 

^ = «M+) + (M_)) /(y) + 2 (M ) %) (179) 

First, let us see if charge particles alone can satisfy Ea. (|178ll . For this to 
be possible we must have dy' f(y') dy" f(y") oc f(y) and f f = 1. 
These conditions are satisfied by 

f{y)^l zSeC ^{y/y) (180) 

where y is a constant specifying the width of the distribution.- 1 This form 
of f(y) is however, too sharply peaked to be consistent with the data. On 
the other hand a Gaussian with the above Oy can approximately satisfy the 
condition dy' f(y') dy" f(y") oc f(y). But in this case, the propor- 
tionality constant at y — is close to 1. All this indicates that we need the 
correlation term involving /o(j/+,y_). 

The condition (|T75|l can be approximately satisfied if /o has the follow- 
ing form 

fo(yi,y 2 ) = q(yi - V2)g((vi + wa)/2) (182) 

where q(y) is a sharply peaked function at y = with a small width A, and 
g(y) ~ (1/ (N))dN I 'dy . Then it can be also shown that h(y) « g(y) and 
1 rv 

dy' dy"f (y',y")*CAh(y) (183) 



•iThis can be easily verified using 

dtanh(x) , o , . , , . 

— = 1 - tanh 2 (x) = sech 2 (x) 

dx 



(181) 
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with some constant C < 1 . In fact the authors of Ref?^ argued that their 
proton-proton result is consistent with having only the neutral clusters. In 
our language, that corresponds to (M±) = 0. 

What the above consideration implies for our charge fluctuation is some- 
what unclear. To the authors' knowledge, direct measurement of charge fluc- 
tuation in the sense defined in the previous sections has not been carried 
out in proton-proton experiments. What have been measured are the C a p 
functions defined in Eq. (|69[1 . In ReJ^I, the maximum height at y 1 = y 2 = 
are given as 

C++ (0,0) w 0.36 
C__ (0,0) w 0.25 

C+_(0,0) « 0.5 (184) 

Using also the overall shape given in the same reference and the facf^that 
(N c h) ~ 6 at this energy one can then infer that 



0.5 ~ 0.6 (185) 



A)/ 



within — 1 < y < 1 without any charge conservation corrections. Correct ing 
for charge conservation is a difficult task to perform. If as asserted in Ref. 
all particles are produced via neutral clusters, then there is no correction 
to perform. If however, some charged particles are emitted independently, 
then this is the lower bound. A more thorough analysis of the available data 
can undoubtedly yield more accurate estimate of the charge fluctuation. 
However, that is clearly beyond the scope of this review. We just make a 
remark here that our estimate of the charge fluctuation in QGP scenario 
(SQ 2 ) Ari / (N c h) An ~ 0-25 ~ 0.3, seems to be still substantially smaller 
than the above proton-proton collision result. 



6.2. Fluctuations in Heavy Ion Reactions 

Unfortunately at the time of the writing of this review, very few published 
data on fluctuations in heavy ion collisions are available. Quite a few pre- 
liminary results are being discussed at conferences, which we will briefly 
mention. But we feel that an in depth discussion of these results prior to 
publication is not appropriate. The pioneering event-by-event studies have 
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been carried out by the NA49 collaboration. They have analyzed the fluc- 
tuations of the mean transverse momentum ^ and the kaon to pion ratioEI. 
Both measurements have been carried out at at the CERN SPS at slightly 
forward rapidities. In Fig. 0] the resulting distributions are shown together 
with that from mixed events (histograms). In both cases the mixed event 
can essentially account for the observed signal, leaving little room for gen- 
uine dynamical fluctuations. Specifically, NA49 gives & Pt = 0.6 ± 1.0MeV k 
for the transverse momentum fluctuations, which is compatible with zero. 
For the kaon to pion ratio they extract a width due to non-statistical fluc- 
tuations of <J n on B tat = 2.8% ± 0.5%, which would be compatible form the 
expectations of resonance decays'"^. 




Fig. 4. Results for the fluctuations of the mean transverse momentum (left) and kaon 
to pion ratio (right). Both results are from the NA49 collaboration. 



The PHENIX collaboration recently also reported their results on mean 
transverse momentum and energy fluctuations at RHIC energies . Similar 
to NA49, their result is compatible with the statistical fluctuations only, 
leaving no room for significant dynamical fluctuations. Their measurement 
was around mid-rapidity with a small azimuthal acceptance of A$ = 58° 

In contrast to that preliminary results by both the CERES collaboration 



k For the definition of Q Vt see previous section. 
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EH at the CERN SPS, as well as the STAR collaboration EH at RHIC re- 
port significant dynamical fluctuations of the mean transverse momentum. 
Both measure at mid-rapidity. It remains to be seen if these difference in 
measurements can be attributed simply to the different acceptance regions 
covered. 

As far as the charge fluctuations are concerned, first results have been 
reported. PHENIxElat RHIC which measures with a small rapidity accep- 
tance, finds c har ge fluctuations consistent with a resonance gas. CERES 
EH and NA49E3 which both measure at SPS energies, report preliminary 
results on charge fluctuations, which are consistent with a pure pion gas. 
However, at the SPS the overall rapidity distribution is rather narrow, so 
that the correlation effect of the resonance gets lost when correcting for 
charge conservation E3 But certainly, non of the measurements is even 
close to the prediction for the QGP. These findings have prompted ideas, 
that p oss ibly a constituent quark plasma, without gluons, has been pro- 
duced However, the measurement of additional observables would be 
needed in order to distinguish this from a hadronic gas. 

But maybe the present range of Ay is so small, that the charge fluctua- 
tions have time to assume the value of the resonance gas. As shown in Fig. 
the results for the parton cascade arrive at the predicted value for the 
QGP only for Ay > 3. None of the present experiments has such a coverage 
yet and thus a detailed analysis of D as a function of Ay is needed, before 
any firm conclusions can be drawn. 

The Balance function measurement at y/s — 130 GeV has been reported 
by STAR collaboration^!. Going from peripheral to central collisions, the 
width of balance function steadily decreases. The trend is what one would 
expect if more of the system is filled with a QGP as the collision becomes 
more central. However, since the reduction is only about 20 % going from 
most perpheral to most central, it is not yet clear whether this signals the 
presence of a QGP or more mundane effect such as the strong flow. 
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7. Conclusions and Outlook 

In this review we discussed several aspects of the emerging field of event-by- 
event fluctuations in heavy ion collisions. We have briefly reviewed the ba- 
sics concepts of fluctuations in a thermal system, and we have discussed the 
possible connection between a idealized thermal system and a "real world" 
heavy ion collision. Furthermore we have introduced the variables used in 
order to present fluctuations and discussed the relation between them. Fi- 
nally we have present the presently published data, which are rather limited. 
But this is soon going to change. 

Crucial for the understanding of fluctuations is the connections to corre- 
lation functions. While the concept of fluctuations might be more intuitive 
and thus might lead us to new insights and approaches, the relevant infor- 
mation is encoded in the correlation functions. For example, consider the 
charge fluctuations, which we have discussed in some detail. It is true that 
like-sign and unlike-sign correlation functions encode all the necessary infor- 
mation on the charge fluctuations. However the physical meanings of those 
encoded information can be made clear only through considering charge 
fluctuations in and out of QGP and its connection to the relevant features 
of the correlation functions. But certainly from the observational point of 
view, the information to be extracted are the correlation functions between 
all possible quanta. Once they are available, the relevant observables can be 
readily obtained by folding the correlation functions with the appropriate 
variables 

On important aspect relevant to heavy ion collisions is the finitcness 
of the system. Thus approximations which are valid for a large system, 
such as the grand-canonical ensemble. Appropriate corrections have to be 
applied to either theory or data in order to extract the properties of the 
bulk matter. While this is not new to heavy ion physics, for the case of 
fluctuations these methods still require more refinement. 

Since the study of fluctuations is a new and developing area of heavy ion 
physics, this review can only be some kind of snapshot at a "hopefully" fast 
developing field. By no means can it be comprehensive. Most of the exper- 
imental data are available, if at all, only in preliminary form. Furthermore 
analysis methods and possible new observable based on fluctuations are be- 
ing developed as we write the review. Therefore our idea was to present the 
basic concepts as well as some of the necessary formalism in a consistent 
and understandable way. To which extent we have succeeded is beyond our 



judgment. 
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